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A search has been performed on 683 pb of data collected by the ALEPH detector at the 
i i Large Electron-Positron (LEP), collider at centre-of-mass energies from 183 to 209 GeV look- 

ing for a Higgs boson decaying into four r leptons via intermediate pseudoscalar a particles, 
^ for a Higgs mass range of 70 to 114 GeV/c 2 and an a mass range of 4 to 12 GeV/c 2 . No excess 

Oh above background is seen and a limit is placed on f 2 = a{ ~ e . 1 ~* z 't h \. x (h — > aa) x (a — > 

t + t - ) 2 in the m h ,m a plane. For m h < 107 GeV/c 2 and m a < 10 GeV/c 2 , £ 2 > 1 can be 
excluded at the 95% confidence level. 

1 Introduction 

Direct searches at LEP2 for the standard model (SM) Higgs boson, h, decaying into bb or t + t 
placed a lower bound of 114 GeV/c 2 on the Higgs mass^. However, fits of the SM to electroweak 
precision data suggest a Higgs with a mass within the kinematic limit of LEP. Additionally, a 
small, non-SM-like excess observed at a Higgs mass of around 100 GeV/c 2 in the bb final state at 
t-h LEP and the fine-tuning needed in the minimal supersymmetric standard model (MSSM) have 

led to the consideration of models, such as the next-to minimal supersymmetric standard model 
(NMSSM) 2 , that feature exotic Higgs boson decays and naturally light pseudoscalar particles, 
a. In these models, new decay channels, such as h — > aa, can dominate over h — )• bb and render 
the Higgs boson "invisible" to conventional searches. In particular, the Higgs can decay into 
four SM particles instead of two, via two intermediate a particles. Several of these possible final 
states, such as h — > 2a — > 46, are already highly constrained by existing analyses; see Ref.^, 
for example. For m a < 2m&, however, the decay a —> t + t~ is expected; this process, with the 
Higgs decaying into t + t~t + t~ for a Higgs mass range of 86 to 114 GeV/c 2 , is not covered by 
existing analyses. To investigate this range, the ALEPH data has been revisited. The present 
analysis is described in detail in Ref.^. 
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2 The ALEPH Detector 

A detailed description of the ALEPH detector can be found in Ref.^ and of its performance in 
Ref.®. High momentum resolution is achieved via a large tracking volume immersed in a 1.5 
T magnetic field. An energy-flow reconstruction algorithm measures total visible energy in the 
event by combining measurements from the tracking sub-detectors and the electromagnetic and 
hadronic calorimeters, and provides a list of reconstructed objects (energy-flow objects) which 



are classified as charged particles (which correspond to charged particle tracks, here called 
tracks), photons, and neutral hadrons. These energy-flow objects are the basic entities used in 
the present analysis. 

During LEP2 the machine operated at centre-of-mass energies from 183 to 209 GeV and 
collected data corresponding to a total integrated luminosity of 683 pb _1 . 

3 Signal and Background Samples 

We revived all steps of the ALEPH analysis framework, including the ability to generate sim- 
ulated samples of standard model background and data. We produced 3000 simulated signal 
events (with h — ^ aa followed by a — > t + t _ ) for each of the three Z decay channels considered 
and for each combination of Higgs boson and pseudoscalar masses in the ranges 70 < < 
114 GeV/c 2 and 4 < m a < 12 GeV/c 2 in steps of 2 GeV/c 2 . For the relevant background pro- 
cesses, our samples were either 10-30 or 300-1000 times larger than the data, depending upon 
the process. 

4 Event Selection 

For the mass range considered, the Higgs is produced approximately at rest, and thus the decay 
h — > 2a — > 4r results in a pair of taus recoiling against another pair of taus. For the a mass 
range considered, the decay products of each 2r system will be observed as a highly-collimated 
jet of charged particles. Due to this high level of collimation, individual identification of taus, 
via standard algorithms, would fail. Instead, we used the fact that each tau decays into either 
one charged particle or three charged particles, and we would thus expect each a jet to contain 
two, four or six tracks. We used the JADE algorithm to form jets with a y cu t chosen to merge 
proto-jets up to a mass of mj e t = 15 GeV/c 2 . 

We considered three possible decays of the Z boson, namely Z — > e + e~ , Z — > and Z 

— > vv, and formulated a set of loose selection criteria (convenient to allow comparison of data 
and simulation at an intermediate stage without compromising the blind nature of the analysis) 
and final selection criteria for each of the two Z decay classes considered: Z — > l + l~ (where / = e 
or n) and Z — > vv. 

4.1 z^i+r 

For the Z — > l + l~ channel, four-fermion background processes are prominent. We used ALEPH 
lepton identification algorithms to mask the two most energetic leptons in the event from the 
list of objects clustered by the JADE jet-finding algorithm. The loose selection consisted of 
the following requirements: Two oppositely-charged, isolated leptons; two jets, well-contained 
within the tracking volume (| cos^| < 0.9 ); and the jets and leptons sufficiently isolated from 
each other (| cos#™ n | < 0.95). The final selection consisted of the following requirements: The 
invariant mass of the lepton pair near the Z mass (80 < mi+i-t^ < 120 GeV/c 2 , where 7 
indicates an isolated photon that may have been radiated from one of the leptons and which is 
added to the di-lepton system if doing so corresponds to an invariant mass closer to the Z mass 
than the di-lepton pair alone); missing energy due to neutrinos from tau decays ($ > 20 GeV); 
jets sufficiently separated (| cos^jj < 0); and a signal-like track multiplicity, i.e., each jet must 
contain either two or four tracks. 

4.2 Z^vv 

The Z —7- vv channel represents a larger branching ratio of the Z than the lepton channel, 
and thus drives the analysis. A major background contribution arises from 77 events. The 



loose selection consisted of the following requirements: Modest missing energy and missing mass 
{$ > 30 GeV and rfi > 20 GeV/c 2 ); exactly two jets, well-contained in the tracking volume, 
with a modest invariant mass cut on the dijet system (| cos#j| < 0.85 and rrijj > 10 GeV/c 2 ); 
requirements on the angle of the missing momentum vector with the beam axis and the total 
visible energy in the event, to reject substantial portions of two-photon-initiated and beam 
background events (\cos6 m i ss \ < 0.9 and E v i s > 0.05 Ecm)] an d modest requirements on the 
most energetic jet (Ej 1 > 25 GeV and containing either two or four tracks). The final selection 
consisted of the following requirements: Less than 5 GeV within 30° of the beam axis, to reject 
events with energy deposits in the forward region of the detector; consistency with the Z boson 
decaying to neutrinos ($ > 60 GeV and rfi > 90 GeV/c 2 ); small aplanarity (< 0.05), consistent 
with two back-to-back, highly collimated jets; and a signal- like track multiplicity, i.e., each jet 
must contain either two or four tracks. 

5 Results 

Based upon these selection criteria, our signal efficiency ranged from ~25% to ~50%, depending 
on Z decay channel, Higgs mass, and a mass. We determined that, for the Z — > l + l~ channel, 
we should expect ~3 signal events versus < 0.2 background events, and for the Z — > vv channel 
our expectation was ~11 signal events versus ~6 background events. 

Systematic uncertainties in our Monte Carlo simulation were estimated to be 5% for all signal 
and 10% for background in the Z — > channel versus 30% for background in the Z — > vv 
channel. We found that the background estimate and the number of events seen in data at the 
loose selection agreed within the systematic and statistical uncertainty for all Z channels. 

For the Z — > l + l~ channels, we observed zero events after applying all selection criteria, 
while for the Z — > vv channel we observed two events. These observations are consistent with 
background. 

We place limits on the cross section times branching ratio of our signal process with respect to 
the SM Higgsstrahlung production cross section, £ 2 = - J^ e + e ^_^+fo) x ~ ^ aa ) x ( a — * t + t~) 2 . 
The limits are based upon event counts in three separate track multiplicity bins (corresponding 
to events with two jets where 1) each jet contains two tracks, 2) each jet contains four tracks, 
or 3) one jet contains two tracks while the other contains four tracks) times each of the three Z 
decay channels considered, resulting in nine categories. The resulting joint probability density 
for the event counts is then used to construct confidence intervals using a generalized version of 
the Feldman-Cousins technique^, which incorporates systematic uncertainties in a frequentist 
way®®. Results are shown, for the 95% confidence level, as a function of rrih (for m a = 10 
GeV/c 2 ) on the left in Fig. [l] and as contours within the rrih,m a plane on the right in Fig. [I] 
Note that our selection criteria do not depend on m/, or m a , and thus our upper limits are fully 
correlated. The observed number of events is consistent with a downward fluctuation of the 
background and, as such, our limits on £ 2 are stronger than expected. 



Also shown on the left in Fig. [T] is the effect of these results upon some possible favored 

scenarios in the NMSSM; see Ref.^, Figures 17 and 21 therein. Our limits highly constrain 
scenarios with tan /3 > 3, while scenarios with tan f3 < 2, where there is a larger branching ratio 
of the Z boson into jets, remain unconstrained. 





Figure 1: Left: Observed and expected 95% CL limit on £ 2 as a function of mj, for m a = 10 GeV/c 2 . Also shown 
are some favored points in the parameter space of the NMSSM. Right: Contours of observed 95% CL limit on £ 2 

in the rrih, m a plane. 



6 Conclusions 



We have performed a search for a Higgs decaying into four taus via Higgstrahlung at LEP2, 
for the process h — > 2a — > 4r and Z — > e + e~, fi + [i~ , or vv, using ALEPH data. We observed 
no excess above background, and for rrih < 107 GeV/c 2 and m a < 10 GeV/c 2 , £ 2 > 1 can be 
excluded at the 95% CL, where £ 2 = -^=z^b x (h — > aa) x (a — > r + r - ) 2 . This analysis 

covers a region of parameter space previously unexplored and further constrains models that 
feature light pseudoscalar Higgs particles and non-standard Higgs decays, such as the NMSSM. 
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